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THE VERTICAL MOTIONS OF MONO- ABUNDANCE SUB-POPULATIONS IN THE MILKY WAY DISK 

Jo Bow 1,2,3 , Hans- Walter Rix 4 , David W. Hogg 4 ' 5 , Timothy C. Beers 6 ' 7 , Young Sun Lee 7 , and Lan Zhang 4 

ABSTRACT 

We present the vertical kinematics of stars in the Milky Way's stellar disk inferred from 
SDSS / SEGUE G-dwarf data, deriving the vertical velocity dispersion, a z , as a function of verti- 
cal height \z\ and Galactocentric radius R for a set of 'mono- abundance' sub-populations of stars with 
very similar elemental abundances [a/Fe] and [Fe/H]. We find that all mono- abundance components 
exhibit nearly isothermal kinematics in \z\, and a slow outward decrease of the vertical velocity dis- 
persion: a z (z,R\ [a/Fe], [Fe/H]) « a z ([a/Fe], [Fe/H]) x exp(— (R — R )/7kpc). The characteristic 
velocity dispersions of these components vary from ~ 15 km s _1 for chemically young, metal-rich 
stars with solar [a/Fe], to > 50 km s _1 for metal poor stars that are strongly [a/Fe] -enhanced, 
and hence presumably very old. The mean <j z gradient (dcr z /dz) away from the mid plane is only 
0.3 ± 0.2 km s _1 kpc -1 . This kinematic simplicity of the mono-abundance components mirrors their 
geometric simplicity; we have recently found their density distribution to be simple exponentials in 
both the z and R directions. We find a continuum of vertical kinetic temperatures (oc a^) as func- 
tion of ([a/Fe], [Fe/H]), which contribute to the total stellar surface mass density approximately as 
T,R (al) oc exp(— <j 2 z ). This and the existence of isothermal mono-abundance populations with in- 
termediate dispersions (30 to 40 km s _1 ) reject the notion of a thin-thick disk dichotomy. This 
continuum of disk components, ranging from old, 'hot', and centrally concentrated ones to younger, 
cooler, and radially extended ones, argues against models where the thicker disk portions arise from 
massive satellite infall or heating; scenarios where either the oldest disk portion was born hot, or 
where internal evolution plays a major role, seem the most viable. In addition, the wide range of 
cr^([a/Fe], [Fe/H]) combined with a constant a z (z) for each abundance bin provides an independent 
check on the precision of the SEGUE -derived abundances: 5[ a /Fe] ~ 0.07 dex and Jpe/H] ~ 0-15 dex. 
The slow radial decline of the vertical dispersion presumably reflects the decrease in disk surface-mass 
density. This measurement constitutes a first step toward a purely dynamical estimate of the mass 
profile of the stellar and gaseous disk in our Galaxy. 

Subject headings: Galaxy: abundances — Galaxy: disk — Galaxy: evolution — Galaxy: formation 
— Galaxy: kinematics and dynamics — Galaxy: structure 



1. INTRODUCTION 

The spatial structure and kinematics of stars in the 
Milky Way's disk appear to be complex, varying dis- 
tinctly with the elemental abundances of different sub- 
populations: more metal-rich populations form thin- 
ner, kinematically cooler sub-components, while metal- 
poor and [a/Fe] -enhanced populations form thicker, 
kinematically hotter disk components (e.g., Fuhrmann 
1998; Chiba & Beers 2000; Feltzing et al. 2003; Lee et al. 
2011b; Bovy et al. 2012a). The elemental abundances of 
stars, in particular [Fe/H] and [a/Fe], encode — albeit in 
complex ways — information about their age and birth 
location within the galaxy (e.g., Schonrich & Binney 
2008). The present chemo-dynamical structure of the 
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Milky Way's stellar disk must reflect the combination 
of the orbits on which stars were formed, and the evo- 
lutionary changes that they underwent since, with the 
abundances being the best practical sub-population 'tag' 
that is preserved throughout the lifetime of a star. As 
such, the chemo-dynamical structure of our Galaxy's disk 
holds unique clues towards understanding its formation 
and evolution. Because a good portion of stars in the 
present-day universe live in galaxies comparable to the 
Milky Way, and because the large majority of Milky Way 
stars live in the disk, dissecting and understanding our 
Galactic disk has broad implications for our understand- 
ing of galaxy formation. 

In a recent paper, we have expanded on earlier work 
studying the Galactic disk as a function of elemental 
abundances, by showing that the geometric structure of 
the disk is relatively simple when viewed as a superposi- 
tion of 'mono- abundance' sub-populations (Bovy et al. 
2012a, Blla hereafter) drawn from SDSS / SEGUE 
(Abazajian et al. 2009; Yanny et al. 2009): the density 
of any set of stars chosen to have a narrow range in 
[a/Fe] and [Fe/H] appears to be well-described by a sin- 
gle exponential in both the radial and vertical directions. 
In particular, Blla found that the (single) vertical scale 
height, h z at a given ([a/Fe], [Fe/H]) varies systemati- 
cally with these two abundance parameters, in the sense 
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that populations, which are metal poor and a-enhanced, 
or a-old, are vertically thicker. It had been established 
that the thicker disk components are a-old (Fuhrmann 
1998; Prochaska et al. 2000; Bensby et al. 2003), but our 
recent work shows that the reverse also holds: a-old sub- 
populations form thicker disk components than a-young 
sub-populations. Furthermore, an analysis of the total 
amount of stellar mass in disk components of different 
thicknesses shows a continuous and monotonic distribu- 
tion of disk thicknesses rather than a simple thin-thick 
dichotomy (Bovy et al. 2012b). 

In this present paper, we follow up on Blla by ex- 
ploring the vertical kinematics, p(v z \z, R) or <j z (z, R), for 
mono-abundance sub-populations, in practice for subsets 
of G-type dwarfs from SDSS / SEGUE in a small range of 
[a/Fe], [Fe/H]). The goals of this analysis are two- fold: 
rst, to see how the distinct and simple spatial struc- 
ture of mono-abundance sub-populations is reflected in 
the vertical motions of these sub-components; secondly, 
to lay the ground work for a dynamical analysis to de- 
termine the gravitational potential near the disk plane. 

Beyond the immediate solar neighborhood (~ 100 — 
200 pc, e.g., Flynn & Fuchs 1994, Nordstrom et al. 
2004), the first determination of the vertical velocity 
dispersion profile was carried out by Kuijken & Gilmore 
(1989), who obtained nearly complete line-of-sight veloc- 
ities for a sample of about 500 K-type dwarfs near the 
South Galactic Pole. For this sample with a fairly broad 
metallicity selection function, they found that the veloc- 
ity dispersion profile, cr z (z), rose from ~ 18 km s _1 near 
the mid-plane to ~ 40 km s _1 at ~ 1.3 kpc above the 
plane. Fuchs et al. (2009) used proper motions of M-type 
stars from SDSS to confirm that the velocity dispersion 
shows a quite strong linear increase away from the disk 
mid-plane, finding an even stronger gradient. Bond et al. 
(2010) used photometric metallicity estimates to isolate 
a low metallicity sample, which showed a considerably 
slower rise of a z (z) away from the plane. 

Much of the rise in cr z (z) found in these analyses is 
presumably attributable to a much higher fraction of 
low- metallicity stars well above the mid plane, stars 
which form a kinematically hotter sub-component. This 
would also explain why the sample with a more re- 
stricted metallicity range (Bond et al. 2010) shows a 
shallower rise towards large \z\. And, just recently 
Liu & van de Ven (2012) confirmed this interpretation, 
using preliminary [a/Fe] values for SDSS / SEGUE G- 
dwarf spectra. They show that the vertical dispersion 
profile becomes approximately isothermal, when splitting 
the sample in a number of [a/Fe] bins. In this paper, 
we expand on these studies, by determining a z (z) for 
sub-populations in a narrow range in ([a/Fe], [Fe/H]), 
which turn out to be exquisitely close to isothermal, i.e., 
a z (z\R, [a/Fe], [Fe/H]) « constant. 

When combined with the spatial structure in the 
radial and vertical direction, the kinematics of these 
mono-abundance populations form a basis for explor- 
ing the viability of various internal or external disk 
evolution and disk heating mechanisms. The verti- 
cal density distribution and kinematics of these mono- 
abundance populations can also provide powerful con- 
straints on the gravitational potential perpendicular 
to the plane for two reasons. First, components 



with a simple spatial and kinematic distribution func- 
tion are more easily modeled, and a density pro- 
file i/*(z|[a/Fe], [Fe/H]) oc exp ( - z/h z ([a/Fe], [Fe/H])), 
<j z {z\ [a/Fe], [Fe/H]) w constant is about as simple 
as it could get. Secondly, all mono- abundance sub- 
components feel the same gravitational potential, and 
hence provide extensive mutual checks on the dynamical 
inferences. 

The outline of this paper is as follows. In § 2 we briefly 
summarize the properties of the data set, and in § 3 we 
describe the approach to estimating a z (z\ [a/Fe], [Fe/H]). 
We then present the results of this analysis in § 4, not 
only a z (z\Ro, [a/Fe], [Fe/H]) but also its dependence on 
Galactocentric radius. We also describe how our ob- 
servations that the mono-abundance components are so 
nearly isothermal in the vertical direction and have such 
a strong dependence of the dispersion on ([a/Fe], [Fe/H]) 
can be used to assess the abundance errors of the SEGUE 
data, independent of all other existing constraints. In § 5 
we put these results into context and discuss their impli- 
cations. We summarize our findings in § 6. 

2. SDSS /SEGUE G-DWARF DATA 

The sample on which we draw in this analysis is iden- 
tical to the one used in Blla but we now also use the 
stars' line-of-sight velocities and proper motions explic- 
itly. For an extensive description of the sample, we refer 
the reader to Blla, Yanny et al. (2009), and Lee et al. 
(2011b), and just summarize some salient points here: 
G-type dwarfs are the most extensive of the systemat- 
ically targeted sub-samples in SEGUE to explore the 
Galactic disk; they are the most luminous tracers whose 
main-sequence lifetime is larger than the expected disk 
age at basically all metallicities. Their rich metal-line 
spectrum affords good velocity determinations (~ 5 to 
10 km s _1 ; Yanny et al. 2009), as well as good abundance 
([a/Fe], [Fe/H]) determinations (5 [a/Fe] ~ 0.1,£ [Fe/H ] ~ 
0.2 dex 8 , Lee et al. 2008a,b; Allende Prieto et al. 2008; 
Schlesinger et al. 2010; Lee et al. 2011a; Smolinski et al. 
2011, which we re-examine in this paper). Given the 
(metallicity-dependent) absolute magnitude of targeted 
G-type dwarfs with 0.48 < (g — r)o < 0.55, ranging from 
4 < M r < 6, the distances to the sample stars range from 
0.6 kpc to nearly 4 kpc. Stars somewhat closer to the disk 
plane are sampled by the lines-of-sight at lower Galac- 
tic latitudes, but the effective minimal distance limit of 
the stars (600 pc) implies that the vertical heights below 
one scale height (\z\ < h z ) of the thinner disk compo- 
nents is not sampled by the data. As in Blla, we use 
a signal-to- noise ratio cut of S/N > 15 rather than the 
S/N > 20 recommended by Lee et al. (2011a); we have 
explicitly checked that our results are the same when 
using S/N > 20. 

The specific analysis in this paper requires the height 
above the plane, the Galactocentric radius, and the ver- 
tical velocity and its error for each star. For a star 

at distance d in the direction (Z,b), i.e., with d = 
d • (cos I cos 6, sin I cos 6, sin b) in the Cartesian Galactic 

coordinate system, we have \z\ = \d ■ e z \ and R = 

\Ro + d — z\. We assume that the Sun's displacement 

8 In this paper, we use 5 to indicate observational uncertainties, 
to avoid confusion with the dispersion symbol a. 



3 



from the mid-plane is 25 pc toward the North Galactic 
Pole (Chen et al. 2001; Juric et al. 2008) and that the 
Sun is located at 8 kpc from the Galactic center (e.g., 
Bovy et al. 2009). We ignore distance uncertainties when 
using the spatial position of stars in the sample in the 
analysis below, as the distance uncertainties are much 
smaller than the spatial gradients in the quantities of 
interest. We do, however, propagate the distance uncer- 
tainties into the velocity uncertainties. 

The motions in the z-direction are a combination of 
the line-of-sight velocities v\ os and the proper motions 
in the 6-direction v z (b, d, v\ os ) = v\ os sin b + d fi^ 
we assume that the Sun's vertical motion is 7.25 km 
s _1 (Schonrich et al. 2010). Given that we expect o z < 
20 km s _1 for the coldest disk components, the v z 
errors, 5 Vz1 matter greatly, both because for stars of 
lower signal-to-noise ratio S Vlos ~ 10 km s _1 , and be- 
cause with typical 8^ b ~ 3.5 mas yr _1 (Munn et al. 
2004) the error contribution from the proper motion 

*u*|/xb = 15 km s" 1 • [d/lkpc] • [5^/(3.5 mas yr -1 )] rises 
linearly with distance, becoming dominant beyond 1 kpc. 
The uncertainties in v z for each star are derived by 
marginalizing over the line-of-sight velocity, proper mo- 
tion and distance errors; the distance errors in turn are 
obtained by marginalizing over the color and apparent- 
magnitude errors that enter the photometric distance 
relation (Eqn. A7 of Ivezic et al. 2008, see also Blla), 
which is assumed to have an intrinsic scatter of 0.1 mag 
in the distance modulus. It is important to note that 
for the present analysis, the detailed understanding of 
the sampling function that was painful and crucial for 
Blla and Bovy et al. (2012b) plays no role here, as long 
as one can assume that the vertical velocities of the tar- 
geted G-type dwarfs at a given Z,6, r, g — r, and [a/Fe], 
[Fe/H] have kinematics identical to their untargeted kin 
with those same non-kinematic parameters. 

We begin with a SEGUE G-dwarf sample that has 
about 28,000 stars with acceptably well-determined mea- 
surements, but here we only use those 23,767 stars that 
fall within well-populated bins in the ([Fe/H], [a/Fe]) 
plane (Blla; a bin is well-populated when it contains 
more than 100 stars). FIG. 1 provides a basic represen- 
tation of the data used in the subsequent analysis: the 
top panels show v z vs. R for a broadly selected a-young 
(left; -0.3 < [Fe/H] < 0.25,0.00 < [a/Fe] < 0.25) and a- 
old (right; -1.5 < [Fe/H] < -0.25,0.25 < [a/Fe] < 0.50) 
sample, respectively. The bottom panels show v z vs. z. 
This Figure immediately illustrates that, as expected, 
the vertical velocity dispersion of the a-old sample is con- 
siderably higher than that of the a-young sample, and 
that there are no strong changes of the velocity disper- 
sion with either distance from the plane or Galactocentric 
radius. 

3. MAPPING THE VERTICAL VELOCITY DISPERSION OF 
G-TYPE DWARFS 

We now describe how we estimate the vertical velocity 
distribution profile for each mono-abundance sub- 
population, characterized by a z (z, R\ [a/Fe], [Fe/H]), 
exploiting the discrete information we have: 
{vz,i,8v z ,i\zi,Ri},i — IjNj f° r eacn abundance bin 
([Fe/H] j, [a/Fe] j). Note that we model the velocity 
dispersion, rather than directly measuring it as the 



second moment of the observed vertical velocities 
(v z ,i — ( v z)) 2 /N in spatial bins. This allows us to 
take the individual observational uncertainties 5 Vz i into 
account and to deal with outliers. 

We presume that the vertical velocity distribution 
p(v z \z,R, [Fe/H], [a/Fe]) can be described by 

p(v z \z,R,5 Vzji ,p) =(l-e c ) • ~~ ~ 

V2ircr z {z,R\p, 5 Vzi ) 

' eXP (2^¥^j) 

+ e c - Pb ac k g r(z,R\p), 

(1) 

where p = /([Fe/H], [a/Fe]), p b ack g r(^, R\p) is a sim- 
ple normalized interloper model (reflecting data anal- 
ysis outliers and halo contamination). Specifically, for 
Pbackgr (z,R\p) we choose a normalized Gaussian with a 
dispersion of 100 km s _1 , convolved with the observa- 
tional uncertainty S Vz i? and e c is a free parameter that is 
the fraction of stars deemed interlopers. In each mono- 
abundance bin, the velocity distribution is presumed 
Gaussian with zero mean and a dispersion a z (z, R\p), 
described by a quadratic function in \z\, as we have no 
strong priors on the functional form of the vertical disper- 
sion profile, except that it be smooth. We assume that 
the radial dependence of the velocity dispersion separates 
from the vertical dependence, and that it is exponential: 

<tI(z,R\p,6 VZi .) =a 2 z (z,R \p u p2,Ps) 

x exp [—2p4-(R — Ro)} 

(2) 

where p^ is the inverse scale length R a _1 for the ra- 
dial change in the vertical velocity dispersion. pi,P2, 
and ps describe the vertical profile of the velocity dis- 
persion, with pi = ^(^1/2), P2 = d& z (zi/2)/dz, and 
p 3 = d<j^(zi/ 2 )/cl^ 2 ; Z1/2 is the median vertical height 
above the plane for each mono-abundance data sam- 
ple. The addition of the observational uncertainty 5% z . 
is due to the convolution of the model with the Gaus- 
sian uncertainty model. Assuming uninformative pri- 
ors on p for each mono-abundance sub-population, we 
then derive, through ensemble-based MCMC sampling 
(Goodman & Weare 2010; Foreman- Mackey et al., 2012, 
in preparation), the posterior probability distribution 
function (PDF) for the parameters of the velocity dis- 
persion p([Fe/H]j, [a/Fe]j) and the contamination frac- 
tion, e c ([Fe/H] j, [a/Fe] j), in each mono-abundance bin 
([Fe/H],-, [a/Fe],). 

To cast these parameterized models in the more intu- 
itive terms of a z (z : R\ [a/Fe], [Fe/H]), one can then use 
samples from the PDF of p to obtain representations of 
plausible dispersion profiles, e.g., cr z (z) at Ro and some 
([Fe/H], [a/Fe]). In each mono- abundance bin the con- 
tamination fraction e c is less than a few percent, and is 
not discussed further. 
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Fig. 1. — Vertical velocities as a function of Galactocentric radius (top panels) and vertical height (bottom panels) for the cn-young sample 
(left panels; -0.3 < [Fe/H] < 0.25,0.00 < [a/Fe] < 0.25 ) and for the a-old (right panels; -1.5 < [Fe/H] < -0.25,0.25 < [a/Fe] < 0.50) 
sample. The median and 16/84 quantiles of the distribution in bins of width ^ 0.4 kpc (as a function of R) and « 150 pc (as a function of 
z) are indicated by the solid line and the shaded regions. 



4. RESULTS 

4.1. Vertical dispersion profile: a z (z,Ro | [Fe/H], [a/Fe]) 

In FIG.s 2-4 we present the results of the our pro- 
cedure for mapping the vertical velocity dispersion of 
the SEGUE G-dwarf sample divided in narrow bins in 
([Fe/H], [a/Fe]). As described above, this procedure is ro- 
bust to outliers and automatically handles the individual 
uncertainties in the observed vertical velocities. FIG. 2 
shows the main result of this paper, namely that the 
vertical velocity dispersion profiles of mono-abundance 
populations, a z (z, Rq | [a/Fe], [Fe/H]), are nearly con- 



stant with height z at some cr^QFe/H], [a/Fe]) — i.e., they 
are isothermal — and that a z ( [Fe/H], [a/Fe]) varies dis- 
tinctly across chemically-selected sub-populations among 
the different abundance bins. The two bottom pan- 
els of FIG. 2 show the vertical dispersion profiles, 
cr z (z, Ro \pbest) f° r eacn abundance bin, where Pbest rep- 
resents the peak of the PDF for the model parameters, 
and the color coding reflects [a/Fe] (left) and [Fe/H] 
(right), respectively. The lines are drawn over the cen- 
tral 95 percent of the |z|-range of stars in each bin, reflect- 
ing the abundance-dependent distance distribution in the 
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Fig. 2. — Best-fit vertical velocity dispersion profile a z (z) and slope of the velocity dispersion as a function of \z\, evaluated at the median 
vertical height f° r each bin, as a function of [Fe/H] and [a/Fe] (top panels) in bins of width 0.1 in [Fe/H] and 0.05 in [a/Fe]. The 

bottom panels show the best-fit velocity dispersion profile (modeled as a 2nd-order polynomial) as a function of height above the plane for 
each bin in the top panels, color-coded for each point in ([Fe/H], [a/Fe]), according to [a/Fe] (left panel) and [Fe/H] (right panel). The 
velocity dispersion in the bottom panels is evaluated over the range in vertical height that contains 95 percent of the data sample for each 
mono- abundance bin; the median height of the observed data is indicated by a dot in the bottom panels. Each ([Fe/H], [a/Fe]) bin contains 
at least 100 stars, and the width of each line in the bottom panels is proportional to the square root of the number of data points. 



sample. Remarkably, cr z is nearly constant within a few 
km s _1 for all individual ([Fe/H], [a/Fe]) bins, while the 
a z values vary from 15 km s _1 to 50 km s _1 among the 
bins. That variation in <j z (z\ / 2 , i^o | Pbest •> [Fe /H] , [a /Fe] ) , 
where Zi/ 2 is the median vertical height above the 
plane for each mono-abundance data sample, is illus- 
trated in the top left panel of FIG. 2, snowing that 
the vertical velocity dispersion increases towards higher 
[a/Fe] and lower [Fe/H] in a pattern that is qualita- 
tively very similar to the corresponding scale-height map 
of mono- abundance components, /^([Fe/H], [a/Fe]) from 
Blla. The top right panel shows the linear slope of the 
a z (z,Ro | [Fe/H], [a/Fe]) profile, dcr z (z)/dz evaluated at 
Zi/ 2 . The slopes are typically at the level of 1 to 3 km 
s _1 kpc -1 with no recognizable dependence on [a/Fe] 
and [Fe/H]. We will return to this remarkable degree of 
isothermality below. 

While FIG. 2 summarizes the best-fit 
a z (z,Ro | [Fe/H], [a/Fe]), FIG. 3 illustrates the un- 



certainties in the inferred parameters p for all mono- 
abundance bins. The PDFs are summarized by the peak 
and 1-sigma uncertainty ellipses of two-dimensional 
projections of the full five-dimensional parameter space 
for each mono-abundance bin. This figure shows that 
the uncertainty on <j z is typically a few km s _1 , such 
that the slight gap apparent in FIG. 2 around a z = 30 
km s _1 is not significant. This Figure also shows 
that there are no significant correlations among the 
parameters, except for a slight correlation between the 
slope and the magnitude of the dispersion and a strong 
correlation between the slope and quadratic term of the 
vertical dependence of the velocity dispersion. The slope 
and quadratic term are anti-correlated, such that the 
overall change in velocity dispersion over the \z\ -range of 
the sample is small in all cases. The parameters of the 
vertical dependence of <j z are not correlated with the 
inverse scale length i?cr _1 ; we will discuss the inferred 
inverse scale length further below. FIG. 4 shows the full 
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Fig. 3. — Parameter uncertainties for the quadratic velocity dispersion profile fits in each mono-abundance bin: the peak and one-sigma 
uncertainty ellipses for each ([Fe/H],[a/Fe]) bin's PDF for the parameters of the vertical velocity dispersion are shown for 2D projections 
of the 5D parameter space. Clockwise from the top, left panel these are: slope vs. velocity dispersion, slope vs. inverse radial dispersion 
scale length, velocity dispersion vs. inverse radial scale length, and slope vs. quadratic term of the vertical velocity profile. There are no 
significant correlations between the parameters, except for a strong anti-correlation between the slope and quadratic term of the vertical 
dependence of a z (z). 



PDF for the slope and quadratic term of the vertical 
dependence of a z for two mono- abundance bins. These 
full PDFs further illustrate the strong anti-correlation 
between the slope and quadratic term. 

By combining the individual mono- abundance PDFs 
for the slope and quadratic term we obtain a joint esti- 
mate for the vertical dependence of the vertical veloc- 
ity dispersion. The peak of the joint PDF occurs at 
da z (z)/dz = 0.7 ±0.5 km s -1 kpc" 1 , d 2 a z (z)/dz 2 = 
—0.12 ± 0.2 km s -1 kpc -2 , with strongly correlated er- 
rorbars. Assuming no quadratic term, the joint estimate 
for the slope is da z (z)/dz = 0.3 ± 0.2 km s -1 kpc -1 . 
Thus, each individual mono- abundance bin is consistent 
with isothermality at the ~ 10 percent (few km s -1 ) level, 
and the mean slope is consistent with zero at the percent 
level (few tenths of a km s -1 ), when compared with the 
range in vertical velocity dispersion present in the full 
disk sample. 

It is hard to judge from FIG. 3 whether the individual 



slope and quadratic term estimates and uncertainties are 
consistent with the joint estimate, in the sense of these 
individual estimates really being noisy estimates of the 
joint estimate. To answer this question, we calculate for 
each individual mono- abundance (slope, quadratic coef- 
ficient) PDF the quantile of the distribution at which the 
joint estimate lies. In practice, this calculation is done 
by binning the individual PDF in two dimensions, sort- 
ing the binned PDF in ascending order (by flattening the 
2D binned PDF into a ID list), and cumulatively sum- 
ming the resulting list and normalizing the total sum to 
one. This list then contains the quantile of the PDF at 
which each bin in the binned PDF lies, and we can find 
the bin and quantile corresponding to the jointly esti- 
mated slope and quadratic coefficient. We find that the 
distribution of quantiles thus calculated is relatively flat 
with no noticeable dependence on [a/Fe], indicating that 
the individual (slope, quadratic coefficient) estimates and 
uncertainties are consistent with being noisy estimates of 
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Fig. 4. — Individual PDFs for the slope and quadratic term of the 
vertical-velocity-dispersion profile for two mono-abundance bins. 
The grayscale is linear and contours contain 68, 95, and 99 percent 
of the distribution and points beyond the last contour are indi- 
vidually shown. For each mono-abundance bin, the slope and 
quadratic term of the vertical dependence of a z (z) are strongly 
anti-correlated, in such a way as to minimize the total change in 
a z with z. The gray histogram in the right inset of each Figure 
shows the PDF for the slope assuming no quadratic term (i.e., just 
fitting a linear profile to the vertical dependence of a z (z))\ in this 
case the slope is tightly constrained to be near zero. 

the joint slope and quadratic coefficient. 

4.2. Radial dependence of the vertical dispersion: 
°~z{R \ [Fe/H], [a/Fe]) 

FIG. 5 presents the dependence of 
<j z (z,R \ [a/Fe], [Fe/H]) on Galactocentric radius, where 
the fit parameter p^ from equation (2) is represented 




-1.0 -0.5 

[Fe/H] 

Fig. 5. — Variation of the vertical velocity dispersion with Galac- 
tocentric radius. Shown is the best-fit inverse radial scale length of 
the vertical velocity dispersion as a function of [Fe/H] and [a/Fe] 
(left panel) and with error bars as a function of metallicity (right 
panel). Points in the right panel are color-coded using [a/Fe] minus 
the median [a/Fe] as a function of [Fe/H]. Symbol size is propor- 
tional to the square root of the number of data points used in the 
point's ([Fe/H], [a/Fe]) bin. The distribution shown on the vertical 
axis in the right panel is the joint probability distribution for the 
(mean) inverse radial velocity-dispersion scale length. The black 
and gray lines show the mean and standard deviation, respectively, 
of this distribution; this mean inverse radial dispersion scale length 
is 0.14±0.01 kpc" 1 (R a w 7.1±0.5 kpc). 

as an inverse scale length, R~ x . The left panel shows 
R~ 1 ( [Fe / H] , [a /Fe] ) . In contrast to a z ( [Fe/H] , [a / Fe] ) , 
there is no discernible systematic dependence of R~ x on 
[Fe/H] and [a/Fe]. While R~ l is not well constrained 
for any individual mono- abundance component, a 
joint analysis yields a PDF for R~ x characterized by 
0.14 ± 0.01 kpc -1 , corresponding to a typical outward 
decrease of cr z (R) ~ exp (— (R — Rq)/(7.1 ± 0.5 kpc)). 
This outward decrease of the velocity dispersion, also 
seen in the luminosity weighted measurements of other 
disk galaxies, presumably reflects the outward decrease 
of the disk surface-mass density and its corresponding 
vertical force. We will briefly discuss implications below, 
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but defer dynamical analyses of the vertical kinematics 
to a separate paper. 

4.3. The vertical kinematics as a test of the abundance 

precision 

The fact that the data sub-sets we dub 
'mono-abundance components' show so nearly 
an isothermal vertical dispersion profile, 
Kda z (z 1/2 )/d 2 ([Fe/H]Ja/Fe])) < 1 km s" 1 kpc" 1 , 
while <j 2 ([Fe/H], [a/Fe]) itself varies by nearly 40 km 
s _1 among the abundance bins, must set a constraint 
on how cleanly our ([Fe/H], [a/Fe]) bins separate abun- 
dances, i.e., it allows a completely independent check 
on the precision (not, accuracy) of the SDSS / SEGUE 
abundance determinations. As the vertical scale heights 
of these sub-populations vary strongly with abundance 
(Blla), stars with high [a/Fe] and low [Fe/H] will 
always increase in relative density towards larger \z\. As 
a consequence, a sub-sample selected nominally by its 
small [a/Fe] and high [Fe/H] values, will always have 
far greater 'contamination' at large \z\ from stars with 
higher [a/Fe] and lower [Fe/H]. As those contaminants 
have a substantially higher <j z , the dispersion profile 
should rise away from the plane if such contamination is 
important. This effect is manifest in samples with poor 
or no abundance preselection (e.g., Kuijken & Gilmore 
1989; Fuchs et al. 2009). 

To turn the observed isothermality of 
a z (z | [Fe/H], [a/Fe]) into a cross-check on 
SDSS /SEGUE abundances (Lee et al. 2008a,b; 
Allende Prieto et al. 2008; Lee et al. 2011a), we 
proceed as follows. We presume that absent any 
abundance errors, i.e., for the perfect mono-abundance 
sub-populations, the dispersion profiles would be per- 
fectly isothermal. This toy model for a z (z) is motivated 
by the fact that the solution for the vertical Jeans 
Equation in cylindrical coordinates with the dominant 
portion of the disk mass at \z\ < 500 pc (i.e., the 
limiting case of a mass sheet, with constant vertical 
force) and an exponential profile for the tracer density 
(as found in Blla), the solution is cr(z) = constant. In 
this context, any force contribution from a vertically- 
extended mass distribution (e.g., the halo), will lead 
to solutions with a z increasing with height, making it 
unlikely that the error-free dispersion profile slope is 
negative (i.e., falling with height). We then assume 
that the observed c^QFe/H], [a/Fe]) (top left panel 
of FIG. 2) and h z ( [a/Fe], [Fe/H]) (Figure 4 in Blla) 
are sufficient approximations to the true distribution 
of these quantities, and we use the normalizations 
from Bovy et al. (2012b). Viewing then the abundance 
uncertainties 5[ a /p e ] an d #[Fe/H] as f ree parameters 
that correlate bins and hence lead to rising dispersion 
profiles, we determine for what abundance errors the 
predicted and observed slopes of the vertical disper- 
sion profile, da z (zi/2)/dz, or shorter a' z , match best. 
Note that error-induced abundance correlations will 
not change the dispersion profile slope in all bins, as 
a z ([a/Fe], [Fe/H]) is a somewhat complex function of 
([Fe/H], [a/Fe]). Specifically, we calculate the PDF for 



'[a/Fe] 



exp 



and 5 [Fe/H] based on an objective function oc 

} ' 5 {°'z,obs - a z,mod(^/Fe]^[Fe/B])) /(<M >oba ) 2 



where <r f Zimod (S[ a /Fe],5[Fs/ir\) is the dispersion profile 
slope predicted by the various assumptions on the errors. 

FIG. 6 shows the abundance-error- induced slopes for a 
few different magnitudes for the abundance uncertainties 
and compares these to the observed slopes of FIG. 2, 
repeated as the top right panel in FIG. 6. It is clear that 
the abundance uncertainties can certainly not be much 
bigger than the SEGUE-reported uncertainties without 
inducing large slopes and clear patterns in the abundance 
dependence of the slope. 

FIG. 7 gives the quantitative results of comparing 
the abundance-errors-induced slopes with the observed 
slopes. The left panel shows the PDF for 5[ a /p e ] an d 
^[Fe/H] using fiat priors on both, while the right panel 
shows the PDF assuming Gaussian priors on the error 
estimates around their nominal value 5[ a /F e ] =0.1 dex 
(Lee et al. 2011a) and 5 [Fe/H] = 0.2 dex (Smolinski et al. 
2011), both with a standard deviation of 0.1 dex. The 
left panel shows that a very clean metallicity separa- 
tion (i.e., very small [Fe/H] errors) are suggested by the 
data. If we include the nominal error priors, we still 
find that the remarkable near-isothermality of the mono- 
abundance dispersion profiles implies that the abundance 
precision are somewhat better than the nominal values, 



with S 



[Fe/H] 



0.15 dex and S\ 



a/Fe] 



0.07 dex. Note, 



that this analysis only speaks to the error precision, i.e., 
the ability to rank stars in their relative abundance, not 
to the abundance accuracy on any absolute scale. 

5. DISCUSSION 

5.1. The simplicity of mono- abundance sub-populations 

The results in this paper show that the vertical motions 
of mono-abundance sub-populations of the Milky Way's 
disk are extremely simple: they are vertically isothermal 
to a remarkable degree — we detect no systematic depar- 
ture from isothermality at the percent level — and expo- 
nentially declining with radius, in a manner that appears 
consistent with the decrease in disk surface-mass density 
(see below). Combined with the findings in Blla that 
the density of these sub-populations is well described by 
a single vertical exponential and a single radial exponen- 
tial, this shows that the phase-space structure of mono- 
abundance sub-populations is about as simple as it could 
be. Chemically-homogeneous sub-components of the 
disk are individually simple disk components that come 
together to create the complex disk structure that is ob- 
served. Typically, this structure has been decomposed — 
in light of this paper and of Bovy et al. (2012a,b)'s re- 
sults, perhaps spuriously — into a small number of sub- 
components with internal metallicity and velocity gradi- 
ents. 

The fact that the mono-abundance sub-populations 
are so close to vertically isothermal is somewhat sur- 
prising, given the inherent complexity in the simple 
([Fe/H], [a/Fe]) abundance plane induced by the (pre- 
sumably) very different star formation histories in dif- 
ferent parts of the disk, and evolutionary processes like 
radial migration, that radially mix stars with very differ- 
ent origins. We expect each pixel in the ([Fe/H], [a/Fe]) 
abundance plane to contain stars with a variety of ages 
and birth locations, with only a rough mean-age increase 
with increasing [a/Fe]. Whether this a priori reason- 
able picture is inaccurate because mono-abundance sub- 



[Fe/H] 



[Fe/H] 




Fig. 6. — Impact of abundance errors on the isothermality of the dispersion profile in mono- abundance bins: Three of the panels shows 
the expected slope of a z {z) for different magnitudes of the SEGUE abundance errors, when assuming that each mono-abundance bin is 
intrinsically isothermal, with velocity dispersions given by the top left panel of FIG. 2 and spatial distributions measured by (Bovy et al. 
2012a, b). The observed slopes are shown in the top right panel. Large abundance errors induce large slopes by 'mixing' populations of 
different a z with proportions that change strongly as a function of \z\. This is especially the case in the transition region between the 
metal-poor, [a/Fe] -enhanced sub-populations and sub-populations with solar abundances. The fact that the observed slopes are < a few 
km s — 1 kpc - 1 and that there is no discernible trend in the observed slopes in the transition region limits the SEGUE abundance errors to 

0.07 dex (see FIG. 7). 



be <5[F e /H] ~ 0-15 dex and 5[ a /F e ] 

populations are dominated by a single coeval population 
even in the light of star-formation complexity and ra- 
dial migration or whether the Milky Way relaxes to this 
simple state, is a question that is best confronted by sim- 
ulations or the next generation of spectroscopic surveys 
that will determine elemental abundances beyond [Fe/H] 
and [a/Fe], such as APOGEE (Eisenstein et al. 2012). 

The velocity dispersions of the mono-abundance sub- 
components increase smoothly from ~ 15 km s _1 for 
the a-youngest populations to > 50 km s _1 for the a- 
oldest. These extrema are typical for thin-thick disk 
decompositions, but the existence of intermediate pop- 
ulations ~ 30 km s _1 is unexpected in the traditional 
picture with a thin-thick disk dichotomy. We apply the 
procedure of Bovy et al. (2012b) for calculating the to- 
tal stellar surface-mass density in each mono-abundance 
bin, correcting the observed number counts for spectro- 



scopic and stellar-population-sampling selection effects 
using stellar-population-synthesis models (we refer the 
reader to that paper for more details on this procedure). 
This stellar mass is not measured dynamically from the 
vertical kinematics. This results in the stellar surface- 
mass- weighted distribution of the of the vertical velocity 
dispersion or kinetic temperature, shown in FIG. 8. To- 
gether with Figure 2 in Bovy et al. (2012b), this shows 
that the distribution of stellar mass in the Milky Way 
is monotonically, and seemingly continuously declining 
from cool, thin disk components to hotter, thicker disk 
components, with no gap or bimodality. This, combined 
with the isothermality of each mono-abundance compo- 
nent, shows that the intermediate populations cannot 
be a mix of the thinner and thicker components, but 
that they are real and contribute an amount of stel- 
lar mass that is intermediate between the thinner and 
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with Gaussian prior : 
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: \ s [a/Fe] = o.i ±o.i h 
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Fig. 7. — Inferred SEGUE abundance uncertainties assuming 
that mono-abundance sub-populations are intrinsically isothermal. 
The left panel shows the PDF for the abundance uncertainties ob- 
tained from comparing the observed slopes with the abundance- 
error induced slopes (see § 4.3 and FIG. 6); small uncertainties are 
preferred. The right panel shows the PDF when applying a broad 
prior centered on the SEGUE-reported abundance uncertainties. 
The preferred abundance uncertainties are still somewhat smaller 
than the SEG UE-reported errors. We infer that <5[ Fe / H ] ~ 0.15 dex 
and <5[ a /Fe] ~ 0.07 dex. 

thicker components as well. This continuity rejects the 
notion of a dichotomy between a thin and a thick disk in 
the Milky Way. 

We can compare our observational results to some sce- 
narios for creating the thicker disk components, in partic- 
ular those where a pre-existing thin disk is dynamically 
heated by an infalling satellite. On this basis, our results 




a\ [km 2 s" 2 ] 

Fig. 8. — The vertical temperature distribution of the stellar 
disk at rsj R ; Distribution of stellar surface-mass density at the 
Solar radius T,ji Q (a^) as a function of vertical velocity dispersion 
squared a\. The thick black histogram shows the total stellar 
surface-mass contributions in bins in cr^, calculated by summing 
the total stellar masses of sub-populations in [Fe/H] and [a/Fe]. 
The stellar surface-mass densities of individual abundance bins in 
[Fe/H] and [a/Fe], calculated by correcting the observed number 
counts of SEGUE G-type dwarfs into stellar masses using stellar- 
population-synthesis models as described in Bovy et al. (2012b) 
are shown as dots, with values for S_R ([Fe/H], [a/Fe]) on the y- 
axis. The points are color-coded by the value of [a/Fe] in each 
bin and the size of the points is proportional to the square root 
of the number of data points that the kinematic fits are based on. 
The temperature distribution shows no evidence for a distinct hot 
("thick") disk component. 

reject massive satellite infall as the dominant disk heat- 
ing mechanism, as such simulations typically display a 
sharp break in the velocity dispersion between the simu- 
lated thin and thick disk (Abadi et al. 2003; Brook et al. 
2004; Villalobos & Helmi 2008) that is in disagreement 
with our FIG.s 2 and 8. This, in addition to the density 
structure in Blla, points toward an internal mechanism 
for creating thicker disk components. Such an internal 
mechanism could be radial migration (Loebman et al. 
2011). The thick disk component could also form in- 
ternally at early times in a turbulent, clumpy disk (e.g., 
Bournaud et al. 2009; Forster Schreiber et al. 2009), al- 
though creating the continuous distribution of vertical 
temperatures requires that the disk remains turbulent 
over a significant fraction of its history (e.g., Forbes et al. 
2011). Whether multiple minor satellite mergers can 
plausibly lead to the continuous trends with elemental 
abundance of the spatial and kinematic structure in the 
Milky Way's disk remains to be seen. 

5.2. The radial dispersion profile and dynamical 
modeling 

In the simplest model of a razor-thin disk and an expo- 
nential tracer density, the solution of the vertical Jeans 
Equation gives the disk surface-mass density in 
terms of the vertical velocity dispersion <j z and the ver- 
tical scale height h z of the tracer population: oc 
a\jh z . Assuming that the tracer's density does not flare, 
this implies that the mass scale length of the disk is equal 
to half of the exponential scale length of the velocity 
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dispersion, and that the dispersion scale length has no 
abundance dependence. We indeed do not observe any 
abundance dependence in the inferred R~ x in FIG. 5 and 
Blla found no evidence for a flare in the tracers' density 
profile. Thus we obtain an estimate of the disk's mass 
scale length of ~ 3.5 kpc from our combined estimated 
of the dispersion scale length of R a = 7.1 ± 0.5 kpc. The 
addition of the dark matter halo means that this is likely 
an overestimate, such that we expect the scale length to 
be < 3.5 kpc. 

Rather than trying to estimate the halo correction 
to the disk scale length as estimated above, we de- 
fer a proper dynamical analysis to a separate paper. 
The simplicity of the mono-abundance sub-populations 
makes them excellent tracers of the local vertical poten- 
tial, because their spatial and kinematic simplicity indi- 
cates that their phase-space distribution function is sim- 
ple. This removes the common difficulty in modeling 
and marginalizing over the phase-space distribution of a 
kinematic tracer sample when doing dynamical inference 
(e.g., Bovy et al. 2010). The large number of indepen- 
dent mono-abundance sub-populations that all trace the 
same underlying gravitational potential will permit for 
extensive cross-checks. The abundance independence of 
the inferred dispersion scale length in this paper is a first 
example of this that indicates the fruitfulness of this ap- 
proach. 

6. SUMMARY 

In this paper we did the following: 

• We modeled the vertical velocity dispersion 
a z (z,R) and its spatial dependence of mono- abundance 
([Fe/H],[a/Fe]) slices of the SEGUE G-dwarf sample. 
We found that all mono-abundance components exhibit 
nearly isothermal kinematics in the vertical direction, 
and a slow outward decrease in the radial direction: 

* z (z,R\[a/Fe], [Fe/H]) « <r z ([a/Fe], [Fe/H]) 

x exp(-(i? - i? )/7kpc) . 

The mean vertical <j z gradient is only 0.3 ± 0.2 km s _1 
kpc -1 , such that the mono- abundance populations are 
consistent with isothermality at the level of one per- 
cent, when compared to the range of velocity dispersions 
present in this sample. 

• We find a smooth variation in characteristic velocity 
dispersion from ~ 15 km s -1 for chemically-young stars 
with solar abundances to > 50 km s -1 for metal-poor, 
[a/Fe] -enhanced stars that are presumably old. The 
mass-weighted distribution of vertical kinetic tempera- 
tures (oc al) is a continuous and monotonically-declining 
function, similar to the mass-weighted distribution of 
scale heights found in Bovy et al. (2012b). This, and 
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in particular the existence of isothermal, intermediate- 
dispersion (~ 30 to 40 km s _1 ) sub-populations, shows 
that the mono-abundance sub-populations cannot be the 
mix of a thin, cool disk component and a thick, hotter 
component. The data therefore reject the notion of a 
thin-thick-disk dichotomy. 

• The remarkable degree of isothermality of the mono- 
abundance components also allows us to independently 
verify the SEGUE abundance precision, as it limits the 
range of abundances actually present within a bin: we 
find that 5[ a /p e ] ~ 0-07 dex an d ^[Fe/H] ~ 0.15 dex. 

• The radial dependence of the vertical velocity disper- 
sion is well described by an exponential. The exponential 
scale length of this radial dispersion profile has no obvi- 
ous abundance dependence, as expected when a massive 
disk contributes significantly to the gravitational poten- 
tial. The exponential decline of a 2 z presumably traces the 
radial decline of the disk's surface-mass density, but the 
quality of the data demands careful dynamical modeling 
beyond the scope of this paper. The simplicity of the ver- 
tical phase-space structure of the mono-abundance sub- 
populations will allow detailed dynamical modeling of 
the Milky Way's potential near the disk plane. 
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